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Abstract. We call an n-tuple Q 1 ,... , Q n of positive definite nxn matrices a- 
conditioned for some a > 1 if the ratio of the largest among the eigenvalues of 
Q i,... , Q n to the smallest among the eigenvalues of Q i,... , Q n does not exceed a. 
An n-tuple is called doubly stochastic if the sum of Qi is the identity matrix and the 
trace of each Qi is 1. We prove that for any fixed a > 1 the mixed discriminant of an 
a-conditioned doubly stochastic n-tuple is n°C>e ~ n . As a corollary, for any a > 1 
fixed in advance, we obtain a polynomial time algorithm approximating the mixed 
discriminant of an a-conditioned n-tuple within a polynomial in n factor. 


1. Introduction and main results 

(1.1) Mixed discriminants. Let Q i,... ,Q n be n x n real symmetric matri¬ 
ces. The function det [tiQi + ... + t n Q n ), where t±, ... ,t n are real variables, is a 
homogeneous polynomial of degree n in t±,... . t Tl and its coefficient 

^TL 

(1.1.1) D (Qi ,... , Q n ) — — - 7 ^— det (tiQi + ... + t n Q n ) 

is called the mixed discriminant of Qi,... , Q n (sometimes, the normalizing factor 
of 1/n! is used). Mixed discriminants were introduced by A.D. Alexandrov in his 
work on mixed volumes [A138], see also [Le93]. They also have some interesting 
combinatorial applications, see Chapter V of [BR97]. 

Mixed discriminants generalize permanents. If the matrices Q i,... , Q n are di¬ 
agonal, so that 


Qi = diag (aa, • • • , «») for i = 1,... , n, 
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then 


(1.1.2) D (Qi, • • • , Q n ) = per A where A = (a^) 
and 

n 

per A = ^ n a ^h) 
i=l 

is the permanent of an n x n matrix A. Here the i-tli row of A is the diagonal of 
Qi and S n is the symmetric group of all n! permutations of the set {1,... , n}. 

(1.2) Doubly stochastic n-tuples. If Q i,... , Q n are positive semidehnite ma¬ 
trices then D (Qi,... , Q n ) > 0, see [Le93]. We say that the n-tuple (Qi,... , Q n ) 
is doubly stochastic if Qi ,... , Q n are positive semidehnite, 

Qi + ... + Qn — I and tr Qi = ... = tr Q n = 1, 

where I is the n x n identity matrix and 1 r Q is the trace of Q. We note that if 
Q i,... , Q n are diagonal then the n-tuple (Qi ,... , Q n ) is doubly stochastic if and 
only if the matrix A in (1.1.2) is doubly stochastic, that is, non-negative and has 
row and column sums 1. 

In [Ba89] Bapat. conjectured what should be the mixed discriminant version 
of the van der Waerden inequality for permanents: if (Q i,... , Q n ) is a doubly 
stochastic n-tuple then 

ri\ 

(1.2.1) D(Qi,... ,Q n ) > — 

n n 

where equality holds if and only if 

Ql • • • Qn I • 

n 

The conjecture was proved by Gurvits [Gu06], see also [Gu08] for a more general 
result with a simpler proof. 

In this paper, we prove that D (Qi,... , Q n ) remains close to n!/n n ~ e~ n if the 
n-tuple (Qi,... , Q n ) is doubly stochastic and well-conditioned. 

(1.3) a-conditioned n-tuples. For a symmetric matrix Q, let A m i n (Q) denote 
the minimum eigenvalue of Q and let A max (Q) denote the maximum eigenvalue of 
Q. We say that a positive definite matrix Q is a-conditioned for some a > 1 if 

^max (Q) — ^Amin(Q) • 

We say that an n-tuple (Qi,... , Q n ) is a-conditioned if 

(1.3.1) A max (Qi) T cxAmin (Q j) foi all 1 ^ j < n. 

In particular, each matrix Qi is a-conditioned, as we allow i — j in (1.3.1). 

The main result of this paper is the following inequality. 
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(1.4) Theorem. Let (Qi, ■ ■ ■ , Q n ) be an a.-conditioned doubly stochastic n-tuple 
of positive definite n x n matrices. Then 

D(Q u ...Q n ) < n^e-C"- 1 ). 

Combining the bound of Theorem 1.4 with (1.2.1), we conclude that for any a > 
1, fixed in advance, the mixed discriminant of an a-conditioned doubly stochastic 
n-tuple is within a polynomial in n factor of e~ n . If we allow a to vary with n then 
as long as a ^jfip- the logarithmic order of the mixed discriminant is captured 
by e~ n . 

The estimate of Theorem 1.4 is unlikely to be precise. It can be considered 
as a (weak) mixed discriminant extension of the Bregman - Mine inequality for 
permanents (we discuss the connection in Section 1.7). 

(1.5) Scaling. We say that an n-tuple (Pi,... ,P n ) of n x n positive definite 
matrices is obtained from an n-tuple (Qi ,... , Q n ) ofnxn positive definite matrices 
by scaling if for some invertible n x n matrix T and real n,... , r n > 0, we have 

(1.5.1) Pi = TiT*QiT for i = l,...,n, 
where T* is the transpose of T. It is easy to check that 

(1.5.2) D (Pi,... , P n ) = (det T) 2 (f[ 7 * j £> (Qi,... , Q n ), 


provided (1.5.1) holds, see [GS02]. 

This notion of scaling extends to n-tuples of positive definite matrices the no¬ 
tion of scaling for positive matrices introduced by Sinkhorn [Si64]. Gurvits and 
Samorodnitsky proved in [GS02] that any n-tuple ofnxn positive definite matri¬ 
ces can be obtained by scaling from a doubly stochastic n-tuple, and, moreover, this 
can be achieved in polynomial time, as it reduces to solving a convex optimization 
problem (the gist of their algorithm is given by Theorem 2.1 below). More gener¬ 
ally, Gurvits and Samorodnitsky discuss when an n-tuple of positive semidefinite 
matrices can be scaled to a doubly stochastic n-tuple. As is discussed in [GS02], 
the inequality (1.2.1), together with the scaling algorithm, the identity (1.5.2) and 
the inequality 

D(Q U ... , Q n ) < 1 

for doubly stochastic n-tuples (Qi,... , Q n ), allow one to estimate within a factor 
of n!/n n ~ e~ n the mixed discriminant of any given n-tuple of n x n positive 
semidefinite matrices in polynomial time. 

In this paper, we prove that if a doubly stochastic n-tuple (Pi,... ,P n ) is ob¬ 
tained from an ct-conditioned n-tuple of positive definite matrices then the n-tuple 
(Pi,... , P n ) is a 4 -conditioned (see Lemma 2.4 below). We also prove the following 
strengthening of Theorem 1.4. 
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(1.6) Theorem. Suppose that (Qi ,... , Q n ) is an a-conditioned n-tuple of n x n 
positive definite matrices and suppose that (Pi,... ,P n ) is a doubly stochastic n- 
tuple of positive definite matrices obtained from (Q \,... , Q n ) by scaling. Then 

D(P 1 ,...,P n ) < n^e-^-V. 

Together with the scaling algorithm of [GS02] and the inequality (1.2.1), The¬ 
orem 1.6 allows us to approximate in polynomial time the mixed discriminant 

4 

D (Q i,... , Q n ) of an a-conditioned n-tuple (Qi, ... , Q n ) within a factor of n“ . 
Note that the value of D (Qi ,... , Q n ) may vary within a factor of a n . 

(1.7) Connections to the Bregman - Mine inequality. The following inequal¬ 
ity for permanents of 0-1 matrices was conjectured by Mine [Mi63] and proved by 
Bregman [Br73], see also [Sc78] for a much simplified proof: if A is an n x n matrix 
with 0-1 entries and row sums 74 ,... ,r n , then 

n 

(1.7.1) per A < (iq!) 1 ^ . 

1=1 

The author learned from A. Samorodnitsky [SaOO] the following restatement of 

(1.7.1) , see also [So03]. Suppose that B = ( bij ) is an n x n stochastic matrix (that 
is, a non-negative matrix with row sums 1) such that 

(1.7.2) 0 < bin < — for all i, j 

Ti 

and some positive integers 74 ,... , r n . Then 

n / .\ 1 In 

(1.7.3) per B < JJ ^-. 

• 1 

i=l 

Indeed, the function B 1 —y per B is linear in each row and hence its maximum 
value on the polyhedron of stochastic matrices satisfying (1.7.2) is attained at a 
vertex of the polyhedron, that is, where bij G (0,1/r*} for all i,j. Multiplying the 
7-tli row of B by n, we obtain a 0-1 matrix A with row sums 74 ,... , r n and hence 

(1.7.3) follows by (1.7.1). 

Suppose now that B is a doubly stochastic matrix whose entries do not exceed 
a/n for some a > 1. Combining (1.7.3) with the van der Waerden lower bound, we 
obtain that 

(1.7.4) per B = e~ n n 0{a \ 

Ideally, we would like to obtain a similar to (1.7.4) estimate for the mixed discrimi¬ 
nants D (Q 1 ,... , Q n ) of doubly stochastic n-tuples of positive semidehnite matrices 
satisfying 

ex 

(1.7.5) A max (Q*) A for i — 1 ,... ,n. 

n 
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In Theorem 1.4 such an estimate is obtained under a stronger assumption that the 
n-tuple (Q i,... , Q n ) in addition to being doubly stochastic is also o-conditioned. 
This of course implies (1.7.5) but it also prohibits Qi from having small (in partic¬ 
ular, 0) eigenvalues. The question whether a similar to Theorem 1.4 bound can be 
proven under the the weaker assumption of (1.7.5) together with the assumption 
that (Q i,... , Q n ) is doubly stochastic remains open. 

In Section 2 we collect various preliminaries and in Section 3 we prove Theorems 
1.4 and 1.6. 


2. Preliminaries 

First, we restate a result of Gurvits and Samorodnitsky [GS02] that is at the 
heart of their algorithm to estimate the mixed discriminant. We state it in the 
particular case of positive definite matrices. 

(2.1) Theorem. Let Q i,... , Q n be n x n positive definite matrices, let H C R n 
be the hyperplane, 

f n 

H = 1 (xi,... ,x n ) : ^Xi = 0 
l i =1 

and let f : H —> R be the function 

f (x 1 ,... ,x n ) = lndet I ^ e Xi Q r 

\i =1 

Then f is strictly convex on H and attains its minimum on H at a unique point 
(£i,... , £ n ). Let S be an n x n, necessarily invertible, matrix such that 



n 

( 2 . 1 . 1 ) s*s = J2 e ^Qi 

i =1 

(such a matrix exists since the matrix in the right hand side of (2.1.1) is positive 
definite). Let 

n = efi 1 for i = 1,... ,n, 

let T — S~ x and let 

Bi = TiT*Qi,T for i = 1,... , n. 

Then (B i,... ,B n ) is a doubly stochastic n-tuple of positive definite matrices. 

We will need the following simple observation regarding matrices Si,... ,B n 
constructed in Theorem 2.1. 
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(2.2) Lemma. Suppose that for the matrices Q i,... , Q n in Theorem 2.1, we have 


tr Qi = n. 

i= 1 

Then, for the matrices B\,... ,B n constructed in Theorem 2.1, we have 

D(Bi,...,B n ) > D (Qi,... ,Q n ). 


Proof. We have 


(2.2.1) D (B 1 ,... , B n ) — (del T) 2 ( Y[n ) D (Q lt ...,Q n ). 

Now, 


q=i 


( 2 . 2 . 2 ) 


and 


II Ti = exp \ & f 1 

i=l v i=1 


-1 


(2.2.3) 

Since (£i,.. 


(del T) 2 — ldet^e €i Qi) = exp {-/ (£i,... , £ n )} 

\ i= l 

, f n ) is the minimum point of / on H, we have 


(2.2.4) /(£i, •••,£») < /(0,... , 0) = lndet Q where Q = yQ i . 

i= 1 

We observe that Q is a positive definite matrix with eigenvalues, say, Ai,... , X n 
such that 

n n 

^ A i — tr Q = tr = n and Ai,... , A n > 0. 

i=l i=1 

Applying the arithmetic - geometric mean inequality, we obtain 

(2.2.5) det <3 = Ai • • • A n < |^ Al + ’^' + An ^ < 1. 

Combining (2.2.1) - (2.2.5), we complete the proof. □ 
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(2.3) From symmetric matrices to quadratic forms. With an n x n sym¬ 
metric matrix Q we associate the quadratic form q : R n —» R defined by 

q(x) — (Qx,x) for x G R n , 

where (•, •) is the standard inner product in R n . We define the eigenvalues, the 
trace, and the determinant of q as those of Q. Consequently, we define the mixed 
discriminant D (g 1? ... , q n ) of quadratic forms qi,... , q n . An n-tuple of positive 
semidefinite quadratic forms qi,... ,q n : R n —* R is doubly stochastic if 

n 

^ qi{x) — ||a;|| 2 for all x G R n and tr qi = ... = trg n = 1. 

i=l 

The property of being a-condil ioned extends to positive definite quadratic forms 
in a natural way. Namely, a positive definite quadratic form if a-conditioned, if 

q(x) < aq(y) for any two x,y G R n such that ||x|| = ||y|| = 1, 

where || • || is the standard Euclidean norm in R n . Similarly, an n-tuple of positive 
definite quadratic forms g l5 ... , q n : R n —> R is a.-conditioned, if each form q ; is 
ct-condit ioned and if 

qi(x) < aqj(x) for all x G R n and all 1 < i,j < n. 

An n-tuple of quadratic forms pi,... ,p n : R n —> R n is obtained from an n- 

tuple qi,... ,q n : R n —> R by scaling if for some invertible linear transformation 

T : R n —> R n and real q,... , r n > 0 we have 

Pi(x) — Tiqi(Tx ) for all x G R n and all i = 1,... ,n. 

One advantage of working with quadratic forms as opposed to matrices is that it is 
particularly easy to define the restriction of a quadratic form onto a subspace. We 
will use the following construction: suppose that qi,... , q n : R n —» R are positive 
definite quadratic forms and let L C R n be an m-dimensional subspace for some 
1 < m < n. Then L inherits Euclidean structure from R n and we can consider the 
restrictions gf,... ,q n : L —» R of gi,... ,q n onto L. Thus we can define the mixed 
discriminant D (gf,... , g m ). Note that by choosing an orthonormal basis in L, we 
can associate m x m symmetric matrices Q i,... , Q m with gi,... , g m . A different 
choice of an orthonormal basis results in the transformation Qi i —> U*QiU for 
some m x m orthogonal matrix U and i = 1,. .. , m, which does not change the 

mixed discriminant D (q i, ... , Q m ) ■ 
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(2.4) Lemma. Let qi, ■ ■ ■ ■ q n '■ IR n —> M be an a-conditioned n-tuple of positive 
definite quadratic forms. Let L C M n be an m-dimensional subspace, where 1 < 
m < n, let T : L —* M n be a linear transformation such that kerT = {0} and let 
Ti,... ,r m > 0 be reals. Let us define quadratic forms pi, ■ ■ ■ ,p m '■ L —> M by 

pfix ) = TiqfiTx ) for x G L and i = 1,... , m. 

Suppose that 

m 

pfix ) = ||x|| 2 for all x G L and \xp l = 1 for i = 1,... , m. 

i=l 

Then the m-tuple of quadratic forms p±,... ,p m is a 4 -conditioned. 

Proof. Since the n-tuple q±, ■ ■ ■ ■ q„ is a-conditioned, we have 

qfix) < aqj(x) for all x G M n and all 1 < i,j < n. 

We define quadratic forms : L —> M, i = 1,... , m, by 

rfix) = qfiTx) for x G L and i = 1,... , m. 

Then 

(2.4.1) rfix) < arj(x ) for all 1 < i,j < m and all x G L. 

Therefore, 

tr rj < a tr r ? for all 1 < i, j < m. 

Since 1 = tr pi = r t tr r,;, we conclude that Ti — 1/tr r* and, therefore, 

(2.4.2) Ti < aTj for all 1 < i, j < m 
Since pi = T,r % , combining (2.4.1) and (2.4.2), we obtain 

(2.4.3) Pi(x) < a 2 pj(x) for all iGf and all 1 < i.j < m. 

Seeking a contradiction, suppose that 

Pj(x) > a 4 pj(y) for some x,y G L such that ||a;|| = ||y|| = 1 
and some 1 < j < m. 

Then, applying (2.4.3) twice, we obtain 

Pi(y) < o 2 p j (y) < a~ 2 pj(x ) < pfix) for all i, 
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so ill the end 


Pi(y) < Pi(x ) for some x,y G L such that ||x|| = ||y|| = 1 
and all i — 1,... , m, 

which is a contradiction since 

m m 

1 = ^Pi(y) = 

i=1 i=l 

This proves that 

Pj{x) < a 4 pj(y) for all x,y G L such that ||x|| = ||y|| = 1 
and all 1 < j < rn 

and hence concludes the proof. □ 

(2.5) Lemma. Let qi, ■ ■ ■ ■ q n '■ R n —* M be positive semidefinite quadratic forms 
and suppose that 

q n {x) = {u,x) 2 , 

where u G R n and ||u|| = 1. Let H — u 1 - be the orthogonal complement to u. Let 
gi,... , q-n-l '■ H —* M be the restrictions of qi,... , q n -\ onto H. Then 

D{qi, ■ ■ ■ , q n ) = D (gi,... , g„_i). 

Proof. Let us choose an orthonormal basis of R n for which u is the last basis vector 
and let Q i,... , Q n be the matrices of the forms gi, ... ,q n in that basis. Then the 
only non-zero entry of Q n is 1 in the lower right corner. Let Q i,... , Q n -1 be the 
upper left (n — 1) x (n — 1) submatrices of Qi,... , Q n - 1 - Then 

det + ... + tnQif) tfi det ^t\Q\ T • • • T tn— i Qn — 

and hence by (1.1.1) we have 

D (Qi ,... , Q n ) — D (Qi ,... , Qn-l j • 

On the other hand, Q i,... , Q n -i are the matrices of gi, ... , q n -i- □ 

Finally, the last lemma before we embark on the proof of Theorems 1.4 and 1.6. 
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(2.6) Lemma. Let q : R n —* R be an a-balanced quadratic form such that tr q = 
1. Let H C R n be a hyperplane and let q be the restriction of q onto H. Then 

a _ 1 

1-< tr q < 1-. 

n an 


Proof. Let 


0 < Ai R ... R A n 


be the eigenvalues of q. Then 

n 

Ai = 1 and A n < oAi, 

i =1 

from which it follows that 

1 a 

Ai > — and X n < —. 
an n 

As is known, the eigenvalues of q interlace the eigenvalues of q, see, for example, 
Section 1.3 of [Tal2], so for the eigenvalues /ri,... ,p n -i of q we have 

Ai < pi < A 2 < ... < A n _i < p n -i A A n . 


Therefore, 



an 


n n— 1 n— 1 

- A* - tvq = J2Pi > - 1_ “- 

i=2 i=1 i =1 


□ 


3. Proof of Theorem 1.4 and Theorem 1.6 

Clearly, Theorem 1.6 implies Theorem 1.4, so it suffices to prove the former. 

(3.1) Proof of Theorem 1.6. As in Section 2.3, we associate quadratic forms 
with matrices. We prove the following statement by induction on m — 1,... ,n. 

Statement: Let qi, ■ ■ ■ , q n ■ —* R be an ct-conditioned n-tuple of quadratic 

forms. Let L C R n be an m-dimensional subspace, 1 < m < n, let T : L —)■ R n be 
a linear transformation such that ker T = {0} and let rj,... , r m > 0 be reals. Let 
us define quadratic forms pi : L —>■ R, i = 1,... , m, by 


Pi(x ) = Tiqi[Tx) for x G L and i = 1,... , m 
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and suppose that 


m 

Pi(x ) = ||a;|| 2 for all x G L and tr pi — 1 for i — 1,... , m. 

1=1 


Then 

(3.1.1) 


D(pu • • • ,p m ) < exp 



i)+a 4 

k =2 



In the case of m = n. we get the desired result. 

The statement holds if m = 1 since in that case D(pi) = det pi = 1. 

Suppose that m > 1. Let L c M" be an m-dimensional subspace and let the 
linear transformation T, numbers r, and the forms p t for i =1.... . rn be as above. 
By Lemma 2.4, the m-tuple pi, ■ ■ ■ ,Pm is a 4 -conditioned. We write the spectral 
decomposition 

m 

Pm{x) = 'y ^ X j (uj , X ) , 

3 = 1 

where u±,... , G L are the unit eigenvectors of p m and Ai,... , A m > 0 are the 
corresponding eigenvalues of p m . Since tr p m = 1, we have Ai + ... + A m = 1. Let 
Lj — ujLj C L. be the orthogonal complement of Uj in L. Let 

Pij : Lj —» R for i — 1,... , m — 1 and j = 1,... , m 


be the restriction of pi onto Lj. 
Using Lemma 2.5, we write 


(3.1.2) 


Let 


D{Ph ■ ■ ■ j Pm) = ^XjD (pi,... iPm-X^Uj^x) 2 ) 

3 = 1 
m 

= X -i D (ui> • • • >P(m-i)i) where 

3 = 1 
m 

y A, = 1 and A j >0 for j = 1,, m. 
i=i 

Oj = trpy + • • • + trp( m _ 1)j - for j = 1,... , m. 


Since 


m — 1 

y Pij{x) — ||a:|| 2 — Pmj{x) for all x G Lj and j = 1,... , m 

i=1 
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and since the form p m j is cr-balanced, by Lemma 2.6, we have 

„4 

(_ 

, . m _ 9 _L _ 

'j 


(3.1.3) 

Let us define 


1 or 

to — 2 -j - - — < (7j < m — 2 H - for j — 1 ,... , to . 

m 


a 4 m 


m — 1 


= 


p,; ? - for i = 1,. .. , m — 1 and j = 1,... , to. 


Then by (3.1.3), 


D (Plji ■ ■ ■ i P{m — l)j) 


CT,- 


to — 1 


m —1 


D (rij ,... , 


(3.1.4) 


< 1- 


+ 


a 


to — 1 to(to 


i \ m ~ 1 

-^J) ) D ( r tr • • • • i.b) 


< expj-l + ^j^^,... ,r (m _ lb ) 


In addition, 
(3.1.5) 


for j = 1 ,... , to . 


tr rij + ... + tr = to — 1 for j = 1,... , to. 


For each j = 1,... ,m, let wij ,... , rf( m -i)j : Lj —)■ M be a doubly stochastic (to — 
l)-tuple of quadratic forms obtained from rij,... , by scaling as described 

in Theorem 2.1. From (3.1.5) and Lemma 2.2, we have 

(3.1.6) D (nj,... ,r( m _i)j) < D (wij,... ,W( m _i)j) for j = l,...,TO. 

Finally, for each j = 1,... ,to, we are going to apply the induction hypothesis 
to the (to — l)-tuple of quadratic forms wij,... : Lj —> M. Since the 

(to — l)-tuple is doubly stochastic, we have 


m — 1 


VJ Wij(x) — ||x|| 2 for all x £ Lj and all j = 


(3.1.7) 


i=l 


and 


tr = 1 for all i = 1,... , to -* 1 and j — 1,... , to. 


Since the (to — l)-tuple wij,... ,is obtained from the (to — l)-tuple 
rij,... by scaling, there are invertible linear operators Sj : Lj —> Lj 

and real numbers //. v ; 7 > 0 for i = 1,... , to — 1 and j = 1.... . rn such that 

Wij(x) — pijrij(Sjx) for all x G Lj 

and all i = 1 ,... , to — 1 and j = 1 ,... , to . 
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In other words, 


, . . _ . Haim— 1)^. . Haim — 1) . _ s 

Wij(x) =HijTij ( SjX ) = —- -Pij ( SjX) = — - -Pi ( SjX) 


°3 




(3.1.8) 


Pij (Tn ~ 1)7* 


Qi ( TSjX ) for all x E Lj 


and all i = 1,... ,m — 1 and j — 1,... ,m. 


Since for each j = 1,... ,m, the linear transformation TSj : Lj —» M n of an 
(m — l)-dimensional subspace Lj C M n has zero kernel, from (3.1.7) and (3.1.8) we 
can apply the induction hypothesis to conclude that 


(3.1.9) 


f m_1 1 1 

D (wij ,... , n; (m _ 1) j) < exp |-(m - 2) + a 4 -j 


for j = 1,... , m 


Combining the inequalities (3.1.2), (3.1.4), (3.1.6) and (3.1.9), we obtain (3.1.1) 
and conclude the induction step. □ 
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